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Abstract 

We present a general formalism to study the growth of dark matter perturbations when dark energy perturbations and 
interactions between dark sectors are present. We show that the dynamical stability on the growth of structure depends on the 
form of coupling between dark sectors. By taking the appropriate coupling which enables the stable growth of structure, we 
find that the effect of the interaction between dark sectors overwhelms that of dark energy perturbation on the growth of dark 
matter perturbation. Due to the influence of the interaction, the growth index can differ from the value without interaction 
by an amount up to the observational sensibility, which provides an opportunity to probe the interaction between dark sectors 
through future observations on the growth of structure. 
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I. INTRODUCTION 



There has been growing observational evidence indicating that our universe has entered the epoch of accelerated 
expansion [1, 2]. Within the framework of Einstein gravity, this acceleration can be attributed to the so called dark 
energy (DE) with negative pressure, which dominates the content of universe at present. The leading interpretation 
of such DE is a cosmological constant with equation of state (EoS) w = —1 [■'>]. Though the cosmological constant is 
consistent with the observation data, at the fundamental level it fails to be convinced: the vacuum energy density falls 
far below the value predicted by any sensible quantum field theory by many orders of magnitude, and it unavoidably 
leads to the coincidence problem, i.e., "why are the vacuum and matter energy densities of precisely the same order 
today?". More sophisticated models have been proposed to replace the cosmological constant by a dynamical dark 
energy in the conjectures relating the DE either to a scalar field called quintessence with w > —1, or to an exotic field 
called phantom with w < —1. It would be fair to say that there is no clear winner in sight to explain the nature of 
DE at the moment. On the other hand, the remarkable discovery of the cosmic acceleration might also be doubted as 
the result of a modification of standard Einstein gravity at large distances. This happens in f{R) theories [4] and in 
braneworld models [^]. How to distinguish such modified gravity theory from that of DE is an imperative task which 
will bring about not only a breakthrough in cosmology, but also in the field of high energy physics. 

The DE is usually not supposed to clump on the scales of the largest cosmic structures, and the most powerful 
way where its nature can be unveiled is to investigate the expansion history of universe[G]. However, the current 
observations of cosmic expansion history cannot break the degeneracies among different approaches trying to explain 
the cosmic acceleration. Considering that the rate of expansion as a function of cosmic time can in turn affect the 
process of structure formation, an alternative channel to study the clustering properties of cosmic structure has been 
used to detect the possible effects of DE. It is expected that the growth history of dark matter (DM) perturbations 
can provide a complementary probe to distinguish the modified gravity from DE[ ]. However, most of these studies 
have neglected the effect of DE perturbations. The DE perturbations do not affect the background evolution, but 
they are crucial in determining the DE clustering properties [!)], which consequently have effects on the evolution of 
DM perturbations. Recently it has been found that for small sound speed and DE EoS prominently away from — 1, 
the presence of DE perturbations may leave a significant imprint on the growth function of DM perturbations [8]. 

In the minimal picture, DM does not feel any significant interactions from DE. Although this picture is consistent 
with current observations, considering that DE accounts for a large fraction of the universe, it is natural, in the 
framework of field theory, to consider its interaction with the remaining fields in the universe. The possibility that 



2 



DE and DM can interact has been widely discussed in [10]-[30]. It has been shown that the couphng between DE and 
DM can provide a mechanism to alleviate the coincidence problem [10, 11, 13, 14]. In addition, it has been argued 
that an appropriate interaction between DE and DM can influence the perturbation dynamics and affect the lowest 
multipoles of the CMB angular power spectrum [15, l(i]. Furthermore it was suggested that the dynamical equilibrium 
of collapsed structures such as clusters would be modified due to the coupling between DE and DM [19, 24]. In the 
presence of such coupling, there has been some concerns about the stability of the perturbations [21]. However, it 
was proved in [22] that the stability of the curvature perturbation depends on the forms of coupling between dark 
sectors. Since observational signatures on the dark sectors' mutual interaction have been found in the probes of the 
cosmic expansion history [20, 25, 26], it is interesting to ask whether the interaction can have an effect on the growth 
of structure, whether it can provide a more consistent check for the coupling between dark sectors. 

In this paper we are going to study the effect of the interaction between DE and DM on the growth function of DM 
perturbations. Recently there is an attempt on this study [27]. Here we will incorporate the DE perturbations in our 
study. We will restrict our investigation to constant sound speed as well as constant EoS of DE. From our formalism 
we will see that the growth of DM perturbations gets more modification due to the dark sectors' interaction than 
that of the DE perturbations. This provides an opportunity to probe the interaction between DE and DM through 
the growth history of DM. 

The organization of the paper is as follows. In the following section we will provide our analytical framework for 
the perturbation equations at the linear level. In Sec. Ill, we will present our numerical results and discuss the growth 
function of DM. We will give our conclusions and discussions in the last section. 

II. ANALYTICAL FORMALISM 

In this section we will derive the second order differential equations for the perturbations of DM and DE with 
couplings between them in a spatially flat Friedmann-Robertson- Walker (FRW) background. The perturbed space- 
time at first order reads. 



ds^ = a^[-(l 4- 2i})dT^ + 2diBdTdx' + (1 + 2(j))5^jdx'dx^ + DijEdx'dx^], 



(1) 



where ip,B,(j>,E is scalar metric perturbations, a is the cosmic scale factor and Dy = {didj — ^Stj'V^). 



We work with general stress-energy tensor 



(2) 
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V.T,';'; = Qf,. (3) 



for a two-component system consisting of DE and DM. Each energy-mementum tensor satisfies the conservation law 

where Q'(^x) denotes the interaction between different components and A denotes either the DM or the DE sector. The 
perturbed energy-momentum tensor reads, 

5v^T(';") = ^{-2[p';, + m{px + pa)]V' + Va + (pa + px)ex + miSpx + Spx) + Hpx + pa)^'} 

= SQl 

d,SV,T>^^^ - ^{[p'x + H{px + Px)]S7'B+[{p'x+p'x)+m{px+px)]ex (4) 
+ ipx + px)V^B' + VHpx + (pa + Px)0'x + (pa + Pa)V2v4 = d.,5Q\x) 

where 6 = V^w, v is the potential of three velocity and the prime denotes the derivative with respect to the conformal 
time r. The perturbed Einstein equations yield, in the linear order 

nV^^-V^(t>' + 2n'^V^B-—V^-B + ]:[V^-fE' = -AnGa^ip + p)e 

a 6 

-d'dj^p - d'dj(j)+^d'djE" +nd'djE' + ^d'djV^E -2Hd'djB - d'djB' = SirGa^Il) (5) 

where Sp is the total energy perturbation, Sp ~ J2\ '^Pa and {p + p)9 = X]a(p^ + P\Wx- 
Considering an infinitesimal transformation on the coordinates [31], 



Sx' = d'p{xn + vl{x^'mvl=0) 



the perturbed quantities behave as. 



(6) 



a 

B = B + 

3 a 

E = E-2P 

V = v + l3' 

9 = 9 + (7) 



Inserting eq(7)into the medium part of eq(4). we obtain the behavior of 6Q^ on the right hand side of eq(4) 

SQo = ,5g" - g"'c" + gV' 

SQp = 6Qp + Q°/3\ (8) 

where 

SQ' = d'SQp + SQlid^SQl - 0). 
and SQp is the potential of three vector 6Q^. This is consistent with the results got from Lie derivatives, 

Cs^Q" = (Sx^g;; - g"(5<, 

<5g- = dQ" - Cs-^Q" , (9) 



which shows that g" is covariant. 

We expand the metric perturbations in Fourier space by using scalar harmonics [•']!], 





a 




= {B^ke-p')Yt^ 
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= {E + 2kP)Y,^;^ 




= (61 + 



(10) 



and the perturbed conservation equations eq( 4) read 



6'^ + - wx)5x = -(1 + wx)kvx - 3(1 + wx)^' + (2V' - 5x)^^ + 



, Spx 
-5px 



(vx + B)' + n{l-3wx){vx + B) 



k Spx , w'^ , ^ , , , a^go 



Sx 



1 + U)A Spx l + Wx 



-{vx + B) + ktp 



vx 



PA PX 

a^SQpx 



il + wx)px {l + wx)px 



(11) 



The interaction Q'j^x) "^^^ decomposed into two parts, 



(tota/) "I" -^A 



(12) 



where Uj^'^^^^^^ is the total four- velocity as defined in[31][21], = h'j^Q'^x) ^ spaeial vector and vanishes in back 
ground F'^ ~ 0[22], h'^" — g^" + C^(^otQ;)C^(iota/) projection operator. The perturbed quantities read 



SQl = SQxU^total) + QxSU^total) 
SQpX = {Q\V(total)+ f\)U (total)! 



(13) 



and Eq( 11) can go back to Eq(20) and Eq(21) in [21], if we neglect the anisotropic stress Tr\. 
Constructing the gauge invariant quantities [31], 

+ 1) 
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Aa 


, P'x^x + B 

= (>x 

P\ k 




Vx 


E' 








+ Ql 




= SQpx-Ql^- 





1 



6 ' 



the perturbed Einstein equations eq( 5) become 



px 



k {H^ ~ = AnGa^ ^ (px + Px) Vx 

A 



(14) 



(15) 



where we have assumed that the pressure perturbation of DE is isotropic = 0. 

Using the gauge invariant quantities eq( 14), we can obtain the linear perturbation equations for DM from eq(ll). 



A' 



Pm k 
PmH 



-kV„, - 3$' + 2* 



Pm U 



Pm k Pr) 



V 



-nVm +k^- ^^Vrn 



pm 



Prn Pn 

Considering the pressure perturbation of DE [21, 22] 



Pd Pd K 



(16) 
(17) 



(18) 



where Cg = |^ [j-f is the sound speed in DE rest frame , = ^ is the adiabatic sound speed , and further noting 
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that 5pd = Spd ^ p'dS 1 we have the gauge invariant form of DE perturbation equations 



A' 



p1Y± 

Pd k 



Pd k Pd K 



v^ + n{i-3w)Vd 



1 + w 



k Pd 

CUAd + ^^)-{C!~Cl)P^^ 
Pd k - 



Pd k 



w 



1 + w 



PdU 

Vd + k-^ 



Pd 



pd 



Pd {l + w)pd 



(19) 



(20) 



Inserting eq.(17) into eq.(16)to ehminate Vm, in the subhorizon approximation k » aH, we obtain the second 
order equation for the DM perturbation 



)a;,+(-a. 



Pni Prn Pn 

Similarly for the DE perturbation we have 



Pn 



■)(7^+^^)-A™(^^)' + (^l^^)'-^^^^-fc2*. (21) 



A. 



Sn'C'.Ad - (^)'Arf + {H{1- 3w) 



-[n + mcl-QwH 



Pd 



1 + w Pd Pd 



Pd 



Pd 



Pd 1 + w Pd 



]A'd-ki 



pd N 



-?,{w"H + wn')Ad - k^{l + w)^ 



Pd 

^[7.(1-3.)- 
Pd 



k'CtAd 
w p'd ^ g^Qgj ^ ^ a^5Q^d' y^ 

1 + w Pd Pd Pd 



(22) 



Changing the conformal time into the cosmic proper time, we can rewrite the above second order equations as 

A^, = 



A™ + 2(i7+-^)A, 



Pm 



Pm Pm 0> Pm 



Pm Pm ^ Prn 



Pni a 



(23) 



Ad+{2H + ZHCl -6wH+^- -^^1 = -3{H + H')C^Ad -(^)A, 

Pd 1 + w pdj Pd 



H{1 - Sw) 



1 + w Pd Pd 



" ' ^'^'^J-3HC! + 3wH-^-9l\Ad^ki'-Sk)-'^C^^Ad-k^l + w)^ 
L Pd ) Pd a 



wH + w{H + H^) 



An 



^'Q"' >Hil-Sw) 



w Pd aQ3\ , 



-(— 



Pd L 1 + w Pd Pd } a Pd 

where the dot denotes the derivative with respect to the proper time t. For the convenience in the foUowing discussions, 
we can further rewrite the second order perturbation equations for DE and DM into dimensionless form 



tPlnA„_ 
dlna? 



1 3 , 

- - -w{i - a, 



din Am 
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'IP 



aQlW H d fa^Ql, 



a dlna \ p 



H'^Pm 



+ - [0„A,„ + (1 - nm)Ad] e 



2 aQl, din Am 1 d f aSQ°' 



H Pm dlna aH dlna \ p. 



m 1 ^ — /nA^, 



k6Ql 



mp. 



pm ^-lnA„ 



-lnA„ 



(24) 
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dlna^ 



dlnAc 
dlna 



+ 



1 3 

2 - 2^"(^ " ^" 



dlnAa 
dlna 



'SC: - 6w 



2aQ2 i« 1 dpd \ dlnAd 
H pd \ + w pd dlna J dlna 



= -3 



1 dH 
H dlna 
kSQl 

aSQ^^' 



1 



aH dlna 



Pd 



+ X (-3C,2 + 3z. - 

1 + wpd dlna H Pd ) [ H pd 



-TTH + (1 + ^) o {^m^™ + (1 - f^m)^4 e 

a'^H'^ I 



1 



3 



dw 



dlna 



1 dH 
H dlna 



1 



(25) 



Hpd y 1 + w Pd dlna Hpd J aH dlna \ pd 

In the subhorizon approximation, from the perturbed Einstein equations eq( 15) we can get the "Poission equation" 

3 







-H' {n^A^ + (i ~ n^)Ad} 



(26) 



where we have used Friedmann equation in the derivation. This equation can be used to build the bridge between 
the matter perturbations to the metric perturbations. 

Since the nature of DE and DM remains unknown, it will not be possible to derive the precise form of the interaction 
between them from the first principle. One has to assume a specific coupling from the outset [11, 29, 30] or determine 
it from phenomcnological requirements [12, 25]. From eq(8)and cq(9), we know that Q"^ is a covariant vector, which 
does not need to depend on the four velocity. For the generality, we can assume the phenomcnological description of 
the interaction between dark sectors in the comoving frame as 



-^{SlPm + S2Pd), 0,0,0 
-— 7((5lPm + S2Pd), 0,0,0 



(27) 



where b\,E2 are small positive dimensionless constants and superindex T is the transpose of the vector. Choosing 
positive sign in the interaction, one can ensure the direction of energy transfer from DE to DM, which is required 
to alleviate the coincidence problem [17, 23] and avoid some unphysical problems such as negative DE density etc 
[21, 25]. In the subhorizon approximation k » aH, 

m 



SQ 



01 



'-{6iSpl^ + S2Sp'd) 



m, 



^Qd ^ -4{Si5pi^ + 52Sp'd) 
Sp'd ^ Spd. 



(28) 



Further noting that the gauge-invariant momentum transfer SQ^^^ refers to the intrinsic momentum transfer between 
dark sectors and as explained in [.SI], such intrinsic momentum transfer is due to the collision of particles from 



different fluids. Such collision can produce acoustics in DM fluid as well as pressure which may resist the squeeze of 
the attraction of gravity and hinder the growth of gravity fluctuations during tightly coupled photon baryon period, 
we therefore set SQp„^ w SQ^^ w [22]. This is a choice of interaction and the result should not heavily depend on such 
setting. 

Employing the above interaction form, we can finally arrive at perturbation equations for dark sectors with constant 
EoS of DE 



dlna^ 



dlna 



1 



dlnAj, 



dlna 



(3(Si + 6— — 

r dlna 



3[exp{ln^) - 1] 



H' 



(52 + 35i (52 + 3(5^ /r + (52(^7 + 1 



5.- 



S2dlnAd Ad 

3 r- expiln- — j 

r dlna Am 



+ (1 - ^m)exp{ln——) 



(29) 



d'^lnAa 
dlna^ 



dlnAa 
dlna 



1 3 n o 



dlnAfi , , 3 

-r^ + i^ + ^)o 

dlna Z 



^^expiln—^) + (1 - ilm) 



a^H^ 



1 + U) 1 + w pd 



dlnAd 
dlna 



3S,rexp{ln^)^^ + 3{^ + l){w - Cl) + 3<5i(:^r + r + r') - 3,5i 
Ad dlna H H 



w — Cl + 5ir{\ — exp{ln—^)) 

Ad 



(1 - 3^) 



where r = Pm/pd and the prime denotes d/dlna. 



III. NUMERICAL RESULTS 



,^2 r'2 

l + W 



il + w + Sir + 62) - 3{dir + 62) 



exp{ln-—) 
Cl P'd 



l + W Pd 



(30) 



In this section we present numerical results of solving the above perturbation equations. We concentrate on the 
behavior of the evolution of the DM perturbation. 

It was argued that the presence of the interaction between DE and DM may give rise to dynamical instabilities in 
the growth of structure [■'>-] . In our general form of the phcnomcnological interaction, when we choose the coupling 
between dark sectors in proportion to the DM energy density by setting 62 = while keeping Si nonzero, we do find 
the consistent fast growth of the fluctuations of DM, as seeing in Fig 1. However, when we choose the dark sectors' 
mutual interaction in proportion to the energy density of DE by taking (5i = and (52 7^ 0, we have the stable DM 
perturbation. This result tells us that the stability in the growth of structure depends on the type of coupling between 
dark sectors, which is consistent with the findings in the curvature perturbation in [22]. 

In the following we focus on the interaction proportional to the energy density of DE to keep the stability of the 
matter density perturbation. Initial conditions arc given at the redshift z ~ 3200 where approximately is the time 
of matter-radiation equality, and we take the adiabatic initial conditions and assume zero initial time derivatives of 
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Figure 1: This figure shows the blow up of the DM perturbation wlien the coupUng between dark sectors is proportional to the 
energy density of DM. Afj = 2.41 x 10~® is the amplitude of curvature perturbations from WMAP five-year results and A„i is 
the density contrast of DM. 

DM and DE perturbations. In our analysis we do not allow w < —1. As to k, we choose its value above O.Ol/iMpc^^ 
so that there is large scale structure data on the matter power spectrum [3, 8]. For the sound speed of DE. Cg , we 
restrict it to be positive and smaller than unity. In Fig 2, we show the evolution of perturbations of DM and DE for 
different values of parameters C'^,w,k and the strength of the coupling 62 between dark sectors. 

Without the interaction between dark sectors, we observed that there is the sensible influence of the DE perturbation 
on the evolution of DM perturbation when the sound speed is tiny enough and the EoS of DE substantially deviates 
from —1. In Fig 2 a, we can see this qualitative behavior that for smaller Cg ,the DE perturbation grows. If we take 
w further away from —1, it will grow more and will influence the DM perturbation. In Fig 2 b, it shows that for 
fixed Cg, DE perturbation grows when w deviates from —1. This is consistent with the result in [S]. However, when 
Cg is not so tiny and w close to —1, in the subhorizon approximation, we observed that the influence of the DE 
perturbation is suppressed. For bigger k, DE influence is even smaller, see Fig 2 c. This result will not change when 
the interaction presents. The influence of the interaction between dark sectors can start to appear in the very recent 
epoch. In Fig 2 d, it shows there in the red circle the small deviation caused by the coupling between dark sectors in 
the matter density perturbation at recent time. 

To see more clearly of the influence of different parameters on the growth history of the DM perturbation, we 
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Figure 2; These figures illustrate the behaviors of dark fluctuations in different cases when varying the effective sound speed, dark 
energy EoS,wave number, and coupling. The solid lines are for the DM perturbation while the dotted lines are for the DE 
perturbation. The solid lines in each panel are largely overlapped, except for that part circled in panel (d) 



introduce the growth index 7 with the definition [8] 

7m = (Inf^m)"^ In ^- 



a dA„- 
\,n da 



(31) 



The growth index is generically not constant which was first emphasized and investigated in terms of cosniological 
parameters in [.34] . The growth index has been argued as a useful way in principle to distinguish the modified gravity 
models from DE models [7, 35]. In Fig 3, we show that the influence of the DE perturbation on the growth index for 
the lack of interaction between dark sectors. The red lines in the figure mark the result without DE perturbation. In 
Fig 3 a, we see that for fixed DE EoS w, when Cg decreases, the growth index with the DE perturbation deviates 
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Figure 3; These figures illustrate the behaviors of dark fluctuations in different cases when varying the effective sound speed, dark 
energy EoS,wave number, and coupling. The solid lines are for the DM perturbation while the dotted lines are for the DE 
perturbation. 



more from the result without DE perturbation. This deviation will be even more prominent when DE EoS is further 
away from —1 as shown in Fig. 3 b. We can see clearly that the difference between the growth index with and without 
DE perturbation can be as big as 0.03. Our numerical result further supports what found in [.<]. In Fig 3 c,d, we 
observed that in the subhorizon approximation the value of k does not influence as much as the parameters w, on 
the result of the growth index. 

In Fig 4, we present our numerical results when we incorporate the interaction between DE and DM. The solid 
lines are for the results with DE perturbation, while the dotted lines are for the results without DE perturbation. 
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Figure 4: The growth index behavior when the interaction between DE and DM presents. Solid lines are for the result with 
DE perturbation, while dotted lines are for the result without DE perturbation. 

It is clearly shown that the growth index got more influenced from the interaction between dark sectors than the 
DE perturbation. If we take the best fitting value of 62 from observations of the expansion history of the universe, 
S2 10~^ [25, 26], its influence on the growth index will even overwhelm that of the DE perturbation. Although 
the enhancement of the growth index due to the DE perturbation and the interaction shown in Fig. 3, 4 is clear, the 
available accuracy from the observations such as DUNE[28] etc. is calculated for the ACDM model and may not be 
true for our interacting cosmology, since the current values of QcOi ^co m^y not typically be equal in the two cases as 
argued in [■>?>]. However, this phenomenon is interesting, as it opens the possibility that the future measurement of 
the growth factor may be helpful to reveal the presence of the interaction between DE and DM. 

IV. CONCLUSIONS AND DISCUSSIONS 

The cosmological observations have provided a firm evidence for significant physics beyond standard models. It 
is clear now that the formation of structure in the universe demands DM and the accelerated expansion of our 
universe requires some kind of DE or a significant infrared modification on Einstein gravity. DE and DM are two 
major components of the cosmic energy budget. It is reasonable to explore possible interactions between them in the 
framework of field theory. 

In this paper we have concentrated on the time evolution of the DM perturbation. We have derived general 
equations to describe the perturbations of DM and DE which incorporates the interactions between them. It was 
argued that the interaction between dark sectors might give rise to dynamical instabilities on the growth of structure 
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[32]. However, we observed that this instabihty depends on the form of the couphng. For example, if we choose the 
interaction in proportion to the energy density of DE, we did observe the stable growth of structure. This result is 
consistent with what is found in the curvature perturbation [22]. 

Besides the interaction between dark sectors, we have also discussed the effects of the nonvanishing DE perturbation 
on the evolution of DM perturbation. Usually in the discussion of the growth of structure, the DE perturbation was 
neglected. For the minimal picture without the coupling between dark sectors, the influence of the DE perturbation 
on the growth function of DM perturbation has been examined in [n]. In our general formalism, by taking the 
appropriate coupling between dark sectors which enables the stable growth of structure, we have found that the 
effect of the interaction between dark sectors overwhelms that of the DE perturbation on the growth function of DM 
perturbation. When the DE EoS w is in the vicinity of —1 somewhere abound the best fitted value at the moment, 
the DE perturbation is suppressed, however, when the interaction presents, the growth index can differ from the value 
without interaction by a big amount up to the observational sensibility. This provides an interesting way to probe 
the interaction between dark sectors through the observations on the growth of structure in large scale [25, 26]. 

It would be of great interest to confront our theoretical work to observations to constrain the interaction between 
dark sectors. However, the data available on the growth of structure are still poor and there is still a long way to go 
before we can talk about precision cosmology in this respect. We will leave our work in this direction in the future. 
On the other hand, it is also very interesting to include the interaction between dark sectors to modify the code in 
studying the N-body cosmological simulations on the structure formation. Work in this direction is in progress. 
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